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ABSTRACT 

Context. The coUisional thick target model (CTTM) of solar hard X-ray (HXR) bursts has become an almost 'Standard Model' of 
flare impulsive phase energy transport and radiation. However, it faces various problems in the light of recent data, particularly the 
high electron beam density and anisotropy it involves. 

Aims. We consider how photon yield per electron can be increased, and hence fast electron beam intensity requirements reduced, by 
local re-acceleration of fast electrons throughout the HXR source itself, after injection. 

Methods. We show parametrically that, if net re-acceleration rates due to e.g. waves or local current sheet electric (6) fields are a 
significant fraction of collisional loss rates, electron lifetimes, and hence the net radiative HXR output per electron can be substantially 
increased over the CTTM values. In this local re-acceleration thick target model (LRTTM) fast electron number requirements and 
anisotropy are thus reduced. One specific possible scenario involving such re-acceleration is discussed, viz, a current sheet cascade 
(CSC) in a randomly stressed magnetic loop. 

Results. Combined MHD and test particle simulations show that local £ fields in CSCs can efficiently accelerate electrons in the 
corona and and re-accelerate them after injection into the chromosphere. In this HXR source scenario, rapid synchronisation and 
variability of impulsive footpoint emissions can still occur since primary electron acceleration is in the high Alfven speed corona with 
fast re-acceleration in chromospheric CSCs. It is also consistent with the energy-dependent time-of-flight delays in HXR features. 
Conclusions. Including electron re-acceleration in the HXR source allows an LRTTM modification of the CTTM in which beam 
density and anisotropy are much reduced, and alleviates theoretical problems with the CTTM, while making it more compatible 
with radio and interplanetary electron numbers. The LRTTM is, however, different in some respects such as spatial distribution of 
atmospheric heating by fast electrons. 
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1 . Basic CTTM Properties and Problems 

Since Ide JageJ (1 19641) and Arnoldy et al. (1968) the Collisional 
Thick Target M odel - CTTM (Brown 1971, 1972; Hudson 1971 
iBrownl 1 19731) - of flare hard X-ray (HXR) sources has be- 
come an almost Standard Model of flare impulsive phase en- 
ergy transport and radiation. It offers a simple and reason- 
ably successful description of several basic features of chro- 
mospheric HXR fla res and even some aspects of the distinct 
coronal HXR flares (iKrucker et al.ll20(j8h . These include predic- 
tion/explanation of: footpoi nt sources ; decreas ing HXR source 
height (e.g. Aschwanden etall 120021: iBrown et al. 2002) and 
source area (Kon tar et al.l 120081) with increasing energy; elec- 
tron time-of-flight e nergy-dependent delays in HXR light curves 
d Aschwandenl2004 . 

However, a number of papers (e.g. iBrown et al.l d 19901) ) 
have reviewed problematic aspects of the standard CTTM 
model and aspec ts of recent data [especially from RHESSI 
- iLin et al.l (|2002[) 1 certainly require modification of the most 
basic CTTM involving a single monol ithic loop. These in - 
clude: the motion of HXR footpoints (iFletcher et al.1 12004 ; 
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the smallness of the alb edo component in HXR spectra 
dKontar & Bro wn 2006; Kasparo va et alj2007 l) compared to that 
expected from the strong downward beaming in the CTTM 
( lBrowniri972l) ; the relative time evolution of the heated soft 
X-ray (SXR) plasma emission measure EM(t) and tempera- 
ture T(t), (e.g. 'Horan 1971!; IStoiser et alj|2008"ah . In addition 
the difference between interplanetary and HXR source elec- 
tron spectral indices is inconsi stent with the CTTM predic- 
tion ( iKrucker et alj l2007l l2009l) . These suggest the need for 
more complex models involving e.g. dynamic filamented struc- 
tures, rather than static monolithic ones, and non-collisional 
effects in electron transport. In terms of theory the main 
CTTM problems are the large fractional instantaneous den- 
sity of the electron beam in the corona and the time inte- 
grated total number of electronsjnjected (e.g. Brown & Melros^ 
ri977l:lBenka & Holmanlll994trBenz &Saint-Hilairell2003 ). The 
beam density problem has worsened as est imates of the beam 
(HXR footpoint) a rea have decreased (e.g. iFletcher & WarrenI 
12003 ). though the 'Kontar eLalJ (1200 8') finding that the HXR 
source area increases rapidly with height may a fleviate this. 
These problems arise from three factors (cf iMacKinnonI 
120061): (a) The high beam intensity demanded by the in- 
efficiency of collisional bremsstrahlung compared with long 
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range Coulomb co Uisional heating of the plasm a. This prob- 
lem is worsened (iMacKinnon & Brownl Il989l) when addi- 
tional energy loss pro cesses are inc l uded such as r eturn cur- 
rent dissipation (e.g. lEmslid 119811; IZharkova & G ordovskvy] 
2005 ), Langmuir w ave generation (e.g. 'Hamilton & PetrosiM 



1987t iMelnik etal.|[T99 9; Kontar 2001; Kontar& Reid 2009]) 
masering (e.g. iMelrose & Dulkiil982 ; MacKinnon et al. 1992), 
electron-whistler interaction (e.g. Stepanov & Tsap 2002), 
Weibel instabiUty (e.g. Karlicky 20091), etc; (b) the TTM injec- 
tion assumption that no acceleration occurs in the radiation re- 
gion so that each injected electron radiates only once and for a 
time no longer than its collisional lifetime tcoif, (c) injection of 
the intense beam is assumed to occur from a tenuous coronal 
accelerator 

The term CTTM is in fact used in two different ways in 
the flare literature. Physically a (collisional) thick target is sim- 
ply one in which the radiating electr ons lose all th eir energy 
(collisonally) irrespective of geometry (lBrownlll97ll) . However, 
the term is often used with refere nce to a particular geometry 
(iBrownll 1 9721 1 1 973t iHudsonll 1 9721) where electrons are injected 
downward into the dense chromospheric target after accelera- 
tion in the tenuous corona. Here we will mainly address this ge- 
ometry though our basic considerations of collisional and non- 
collisional transport are relevant to at least some of th e types of 
coronal HXR source reviewed by lKrucker et al.l (1200 8'). Our dis- 
cussion also applies both to HXR footpoi nts of static mo nolithic 
loops and to the scenario described by Fletcher et al.l (l2004i) . 
where footpoint HXR sources move, probably as a result of pro- 
gressive magnetic field line reconnection. 

In this paper we discuss problems (a) - (c) above and pro- 
pose a modified thick target scenario involving similar geometry 
and injection but replacing assumption (b) by a re-acceleration 
process acting within the HXR radiating volume. The resulting 
increase in electron lifetime to much greater than tcoii increases 
the photon yield per electron and reduces the necessary electron 
replenishment rate, beam density, and anisotropy. 



2. HXR Source Requirements 

2.1. Model-independent Nonthermal Emission Measure 

The instantaneous bremsstrahlung output 7(e) (photons s ' per 
unit photon energy e) from a source volume V, with local plasma 
density np(r), and fast electron flux spectrum F(E, r) at position 
r and bremsstrahlung cross section Qsie, E) dififerential in e is 
(Brown 1971) 



7(e) 



F(E)QB(e,E)dE 



(1) 



where the source means are n - J^n(r)dV/V; F(E) — 
J^n(r)F(E,r)dV/(nV). For prescribed Qb, 7(e) is thus related 

unique ly to the quantity nV FjE) regardless of how F{E) is pro- 
duced (iBrown et al.ll2003l) . The spectral shape of 7(e) is fixed 
by the shape of F(E) while the absolute scale of 7 is fixed by 
a spectrum dependent factor of order unity times nVFi where 
Fi - L F(E) is the total mean electron flux above some refer- 

ence energy E - Ei. Following iBrown et al. I (120031) . nVF{E) 
is often used as the fundamental unknown 'source' function 
in inferenc e of HXR e lectron spectra from data on 7(e) ( e.g. 
iPiana etal.l 12003; Ko ntar et al.1 [20041: iMassone et all 120041: ?). 
An equivalent HXR source property which is more readily en- 



visaged physically than the total nVFi is the total nonthermal 
emission measure of electrons of E > Ey given by 



Jv 



ir)n{r)dV - nn\V - n\nV 



(2) 



where ni(r) = jj' F(E,r)dE /v(E) is the local density of elec- 
trons of £■ > El. EMi can readily be used for example to 
find the mean fractional density fi of fast electrons if the to- 
tal (thermal) emission measure EM ^ n V in V is known, viz 
ni/n = EMi/EM. Note also that we can write EMi - nNi 
where Ni is of order the total number of fast electrons in V. 
Numerically, in a typical large event, the necessary EMi is 
> 10'*'' cm""* for El - 2Q keV so any model of an intense HXR 
source must involve conditions satisfying 



EMi = nniV =/in2y 



lO^'^fml^Vn > lO'^^cm-^ 



(3) 



with n = 10'"nio cm-^ V = lO^^VavV cm^ etc. This 
shows that coronal sources alone can only generate large HXR 
bursts if they have unusually large volume and/or density (e.g. 
Veronig & Brown 2004; Krucker et al. 2008). Maintenance of 
this EMi in the CTTM case requires that electrons be injecte d 
at a rate [cf Equation ([7]) below and dBrown &Emsliel[l988l) 1 
Ti > 10^^ s"' above 20 keV. This large value is the origin of : (i) 
the problematically large number of total electrons processed by 
the accelerator during event duration t„, viz. ^ fiTo or around 
10^^ in a few 100 s, equal to lOOx the total electrons in a loop 
of y = lO^^cm-^, n = 10''^ cm"-'; (ii) the high beam density 
«! - TijAvi over area A. For A -2x lO'*" cm^ (or ^2" square) 
this gives «i ^ lO'^cm""*, a density as high as the coronal loop 
plasma density in which the intense CTTM beam propagates. 



2.2. Electron Lifetime and Model-dependent Replenishment 
Rate 

The instantaneous values of nVF{E), EMi etc in practice 
change as the electrons evolve. In cases where electron lifetimes 
T are short (compared to event duration or observational integra- 
tion times) it is necessary to sustain nVF(E), EMi etc and hence 
7, by maintaining the numbers of the electrons of £ > Ei at a 
rate given roughly by 



Ti ^ NilT^EMilinr) 



(4) 



This can be either by injection of fresh electrons from outside 
the HXR source to replace decaying ones (as in tbe CTTM), or 
by a local reaccelation process acting on those inside the source 
to offset their energy losses. The latter option has received very 
little attention in the HXR source literature and is the one we 
focus on in this paper In the case of the CTTM model, main- 
tenance of Ni is by replenishing injection from the corona and 
r - TcTTM here is the electron collision time 
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(5) 



Since tcoii ^ 1/n, by Equation (HJi the injection rate Ti re- 
quired to sustain EMi,F{E) is independent of n (lBrownlll97lh . 
If T is reduced below tcu by non-coUisonal losses then the nec- 
essary Ti is increased and the problems of the CTTM worsened. 
Of much greater interest are situations where the lifetime t in- 
side the HXR source is somehow enhanced over tcu because 
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one can then (Equation dU) attain the same EM\ for a smaller 
replenishment rate Tx but the same instantaneous total N\ . 

While increasing r reduces Tx the consequences for fast 
electron density n\ in the source depend on the geometry of 
their propagation. For example, if the fast electrons were being 
injected into a HXR source from above, increasing r while con- 
taining them in the same V, n (eg by scattering or magnetic trap- 
ping), would leave the fast electron number density unchanged 
but they would last longer and sustain EM\ for smaller T\ . If, on 
the other hand, they propagated freely downward, their longer t 
would cause them to penetrate more deeply, increasing the HXR 
source V and sustaining EM\ with a smaller n\ but larger V. We 
discuss the latter situation again in Section 4. 

The above estimates of the necessary Tx etc in terms of a 
single T value are only approximate. To get a more accurate pic- 
ture of how electron supply requirements are modified by non- 
collisional energy losses and gains it is necessary to look more 
closely at actual photon yield and its relation to electron trajec- 
tories E{t). 

2.3. Electron Trajectories and Photon Yield 

In general the number ^(e) of photons per unit e emitted during 
the lifetime of an electron of initial energy £, is 



((e,Ed 



Jl(E>6) 



n(r(t)Mt)QB(e,E)dt 



where n(r(t)) is the plasma density along the electron path, and 
v(f) = (2£'(f)/'We)'^^ the electron speed while t(E > e) is the total 
of all intervals during which E{t) > e. As electrons tend to decay 
to £ < e (or escape) they have to be maintained at a spectral 
rate ;F*(£',) (s"' per unit £,) to sustain the value of nyF(£') and 
hence 7(e). 7(e), T^jEt) and nV F{E) are inter-related by lBrowrJ 
(1 1 97 Ih : iBrown & Emslid (1 1 988h 



7(e) = X"'r.(£.)^(e,£,)fl'£. 
nV l°°F(E)QB{e,E)dE 



(7) 



where ^(e, £,) is now the mean value for a large number of elec- 
trons of the same initial E - E^ since in general E{t) can differ 
greatly between electrons of the same especially in the case 
of stochastic acceleration - see below. [Even for purely Coulomb 
collisional transport there is dispersion in E{t) for given Et due 
to the spread in impact parameters and the finite thermal speed of 
target particles. Bot h of these are s mall in the CTTM and are usu- 
ally neglected - ('e.g.l Brownll97lt) 1 Clearly the ;F*(£'.) necessary 
for given 7(e) is related to 1 or, crudely, to 1 /t as discussed in 
Section 2.2. Note that when varies on timescales shorter than 
the electron time of flight in the HXR source. Equation (|7]i has to 
be modified to allow for energy dependent time delays between 
features in !F, and in 7 -i.e. acceleration and propagation effects 
are conv oluted in time . This has been discussed in the collisional 
case by Emslid Jl 9831) : lAschwan den ( 2004 ) . It is even more rel- 
evant to the situations discussed here which specifically involve 
extended electron lifetimes t. 

In the CTTM, with radiation only in the collisional propaga- 
tion region (and no acceleration), the mean dEfdt - E - Ec„u - 
-Knv/E where K = 2ne^K with A the Coulomb logarithm. The 
mean E(t) is thus monotonic so the maximum E = E^is the ini- 



tial/injection energy and we can write dt - dEji-E) - dEI\E\ 
and replace the t integration (|6]l by the E integration 



CcTTMie, £.) = i J EQeie, E)dE 



(8) 



It is the small value of eQsie, E) compared with KIE^ here 
that makes collisional bremsstrahlung an inefficient source of 
HXRs in any model, and demands large electron injection rates 
and beam power. Even if we can increase the electron lifetime 
and reduce the necessary number supply rate the power required 
is unchanged or may even be increased. 

Any non-collisional transport process which acts solely to 
add energy losses E to the collisional ones can only reduce ( 
below CcTT M and so increase the necessary and power re- 
quirements [Mac^nno^ (|2006|) . The only processes capable of 
allowing ^ > (cttm, hence reducing are ones which tend on 
average to increase the mean electron lifetimes over tc„ii. (We see 
below that the actual effect of this on ^(e) depends on the form 
of E{t) and of Qb{£, E)). Physically this corresponds to accelera- 
tion inside the HXR source, a process rather arbitrarily excluded 
in conventional CTTM assumptions. The effect on ^ of chang- 
ing E is not immediately obvious as we show by considering 
some simple parametric forms (p{E) to describe the effect of the 
acceleration relative to collisions, viz. 



(6) E^E, 



noncoll "r ^coU 



cf>(E)E, 



coll 



(9) 



To measure the effect of varying on ^ we have to adopt 
a specific form for QB(e,E) and we first consi der the Kramer s 
form Qbk - Qo/^E (e < E) with Q„ a constant (lKramersll923h . 
for which a measure of ( is the quantity 



dE 



(10) 



This simplification lets us give several illustrative analytic 
examples of the dependen ce of ( on trajectories E(t) (cf 
iBrown & MacKinnonlll985l) . The true Qb behaves in a more 
complex way the consequences of which we mention below. For 
collisions only (CTTM), (p - \ and ^ - E,„ax - e. Other informa- 
tive cases are - 

- (i) <p{E) - VE ^ — » oo since the electron formally has 
infinite lifetime. Physically this contrived idealisation would 
be like dragging an electron at constant speed through the 
plasma, energy supply exactly offsetting losses and making 

T — > OO 

- (ii) (piE) - constant C 

(a) C > (net energy loss) ^ ^ = ^ = ^ so that ( 
is only enhanced in this case for < C < 1 which is also 
unrealistic corresponding to to fine tuning of Ea to partially 
offsett losses Ecoii but not reverse them to a net gain. 

(b) C < (net energy gain). Here E(t) increases indefinitely 



{Efj 



oo), as f —> OO and ^ = (-EmnA- - f)/|C| — > oo. Though 



an infinite lifetime is clearly unphysical, arbitrarily increased 
^ is possible if arbitrarily high E,„ax is reached, 
(iii) (f>{E) - -(piiE/Ei)" with (f>\ > (net energy gain). Here 
again there is a formally infinite lifetime with Ei„ax ^ oo as 
r — > oo but, for a \, 



(a - 1)01 



e \ I t-^n 



(11) 
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Fig. 1. Cross section QB(e,E) at e = 20 keV showing how the 
emission contribution per unit e peaks at different E for different 
Ob- (Dotted = Kramers, dashed = Bethe Heitler, solid = exact 
(lHauglll997l) ^ 



This diverges for a < I but is finite Va > 1 despite the infi- 
nite lifetime (E„ax °°)- This is because Qbk ^/E, with 
maximum value at £ = e so that the contribution to ^ falls 
as E increases and the total is finite for any sufficiently fast 
acceleration (a > 1). 

These examples show how f can depend on the specific form 
of the electron trajectory E{t). In addition, ^ depends on the form 
of Qb in relation to E(t). For any e, f will be largest when E(t) 
maximises the time spent near the value of E where Qb peaks. 
For the Kramers Qbk used above this is at £ = e but even for the 
next simplest approximation - the non-relativistic Bethe Heitler 
form Qbbh - the peak is substantially shifted to E ^ 1.7e as is 
also the case for the full cross section as given bv lHaugl(ll997h 
- see Figure 1. Different forms of E(t) convolved with these 
Qb can result in substantial differences in rates of emission and 
hence in ^. 

These special cases illustrate how ^ can be enhanced by re- 
acceleration Ea either by making the net loss rate l^l < \Ecoii\ or 
by creating a net gain rate so that > both increasing the 
electron lifetime at £ > e. This reduces the necessary injected 
beam density and total numbers of electrons by prolonged re- 
acceleration of them in the HXR source after injection. We call 
this the Local re-acceleration Thick Target Model (LRTTM). In 
Section 3 we discuss a specific physical energy release scenario 
where these LRTTM requirements may be met. We emphasize 
again that reducing 'Fi in this way does not reduce the power 
that has to be delivered. This is always at least the value in the 
CTTM since, for every erg of collisional bremsstrahlung emit- 
ted, of order 10^ erg go into long range collisional energy losses. 
However, in the LRTTM, most of that power is delivered in the 
HXR source rather than in an external accelerator/injector of 
electrons as in the CTTM case, a point to which we return in 
Section 4. 

To see whether and to what extent this happens in any par- 
ticular LR scenario we have to recognise that the actual photon 
yield ((e) during the lifetime of an electron in such scenarios is 
more complicated than discussed above. The f integral in equa- 
tion (|6]l cannot be written simply as an integral over E, as it can 
in these cases, since: 

1 . E(t) is no longer monotonic in general, with E taking values 
> 0, < 0, or at different parts of its path. Then the t integral 



can only be written as a sum of E integrals with one for each 
f segment in which E(t) is monotonic (with £ > or £ < 0) 
plus integrals over t itself when £ = so that dt does not 
transform to a finite dE. In practice one reverts to the basic t 
integration (|6]l. 

Even if the change in variable from f to £ is useful, the upper 
limit in the £ integrals is no longer the initial energy £« (as 
it is in the CTTM) but the maximum value E^axiEt) reached 
during the electron lifetime at £ > e 

In re-acceleration, e.g. by waves, the trajectories £(f) are not 
only non-mono tonic but may well be highly stochastic, dif- 
fering between electrons of the same initial £, (cf Section 
3 for a specific example). There is then no well defined de- 
terministic yield ((e, £,) for electrons of initial E - E, and 
the total yield has to be found numerically by evaluating ex- 
pression (|6]l for each elec tron and summing thern , or using 
statistical techniques (e.g. lBian & Browning|l2008b . 



3. Current Sheet Cascades (CSCs) as One Possible 
LRTTM Scenario 

The LRTTM idea that local reacceleration of electrons inside 
the thick target HXR source can greatly increase their photon 
yield by prolonging their lifetimes to » tcoii is a quite general 
one which might be realized for many different (re)acceleration 
mechanisms. The basic requirement is some source of strong 
electric fields distributed through the so urce and this might b e 
achievable in a va r iety o f ways - e.g. iLionello et all j l998l) . 
iFletcher & Hudsoiil ( 1200 8'). In this Section we focus on one pos- 
sibility to illustrate the idea in some detail. 

3.1. Energy Release and Electron Acceleration In CSCs 

The CTTM idea of separation of the acceleration and radia- 
tion volumes had its origins partly in t he idea s that : energy 
is most easily stored in the corona dSweetl Il958h : accel- 
eration is more efficient in a tenuous collisionless volume 
(e.g. Hamilton & Petrosian 1992; Miller etal. 1997), while 
bremsstrahlung gives most volumetric yield at high densities. In 
such cases, magnetic energy release is assumed to be driven by 
organized and continuous twist or shear of larg e scale magnetic 
structures (isolated loops or arcades) - e.g. iForbes & Priest! 
(1995). An alternative is distributed sm all scale release of en- 
ergy in a Current Shee t Cascade (CSC) dGalsgaard & Nordlundl 
119971; lGalsgaardll2002l) . resulting fr om the 3-D MHD response 
of a loop to a random underlying photospheric driver. This 
gives a specific physics-motivated example of the type of local 
re-acceleration scenario discussed schematically in Section 2. 
After a few Alfven times (sees), Lorentz forces create stresses 
along the entire loop and form a hierarchy of reconnecting 
current sheets, leading to plasma jets. These perturb the 
neighboring plasma and eventually create a turbulent current 
sheet cascade (CSC) throughout the volume from large scale 
current sheets (CSs) to numerous small scale CSs in which 
energy is dissipated randomly everywhere. CSs appear and 
disappear over short times but with an overall quasi-steady 
turbulent state. This energy release and acceleration process 
operates not only in the corona but also in the chromosphere e.g. 
(iDaughton et al.ll2008l) . Several papers have discussed particle 
acceleration in the CS C electric fields & in such a dynamic 
enviro nment dAnastasiadis & Vl ahosI 119941; IPmitruk et all 
200 31: lArzrier & Vlahos 2004; Iviahos & GeorgouiisI [2004 
,Vlahos et all 12004 .Turkmani et all I2005L 120061) . the last two 






Fig. 2. Panel a. Examples of trajectories E{E,,t) for nine electrons re-accelerated in the chromosphere, and one (dotted curve) 
undergoing collisional losses only, after arriving from the corona with initial energy = 50 keV. Panel b. Relative photon 
production rate at 20 keV. Panel c. Relative cumulative photon production ^(e, t)dt keV ' at e = 20 keV for electrons shown in 
Panel a. Time is in units of the collision time for a 50 keV electron {tcoii - 0.0032) s. The emission rates in b. and total emissions in 
c. have been divided by the total yield (cttm keV"' at 20 keV of a purely collisional electron starting at = SOkeV. In c, the ((t) 
curves that stop at some t have attained their asymptotic values there, the electron E(t) having dropped below 20 keV thereafter 



of which simulated particle acceleration in the £ fields present 
in the coronal case of such models. For the present paper 
we conducted similar calculations for already energetic elec- 
trons injected into the dense chromospheric part of the loop 
. Trajectories E(t) of test particles, with prescribed initial 
energy £, randomly injected in space, were traced in a frozen 
"snapshot" of the MHD fields since acceleration times (<s; 1 
sec) are much shorter than the overall MHD evolution timescale 
(~ seconds). Electrons gain or lose energy stochastically as they 
travel along guide fields and pass through the local CS £-fields 



and lose it collisionally, following complicated trajectories as 
they receive kicks of various signs and strengths. Some test elec- 
trons encounter few or no CSs and decay purely collisionally, 
rejoining the local thermal plasma. 

Following fTurkmani et alj (|2005LI2006|) we found that, in the 
corona, strong acceleration of a substantial fraction of thermal 
electrons occurs for reasonable values of the resistivity rj so long 
as super-Dreicer fi occur in some of the current sheets. In com- 
mon with most acceleration modelling, it is hard to assess realis- 
tically how high this fraction is in a test particle approach involv- 



6 



Brown et al: Thick target model with local reacceleration 



ing scaled numerical resistivity and resolution. If the fraction 
becomes very high (as it has to in the CTTM to create enough 
HXRs) the validity of the MHD/test particle approach becomes 
questionable since the associated currents should be allowed to 
feed back on the MHD field equations. As we show below, in 
the LRTTM, the necessary 'F is reduced, which alleviates this 
issue. Our goal here is simply to show the implications for HXR 
source requirements if extensive re-acceleration does occur 

After undergoing acceleration in the corona, electrons 
mainly precipitate into the chromosphere where many of them, 
instead of rapid collisional decay, undergo re-acceleration in the 
chromospheric CS & fields. In this paper we therefore only dis- 
cuss what happens to these electrons once injected into the chro- 
mosphere, namely a substantial fraction of them survives at high 
energies for many collision times, increasing the photon yield ^ 
over the purely collisional CTTM value. This is a good example 
of the type of LR scenario suggested schematically in Section 
2 since the CTTM distinction between acceleration and radia- 
tion regions disappears, and (re-)acceleration occurs in the HXR 
source. 

3.2. CSC Simulations of Electron (Re)-Acceleration 

For the coronal part of the loop we adopted plasma values 
n = 10'° cm ^ T = lO^K, B = lO^G and for the chromospheric 
part n = lO'-^cm"^, T = 10'*K, B - lO^G, though we recognise 
that these vary in space (especially n in the chromosphere). The 
depth of the chromosphere is taken to be 7.5 x 10^ km. Here we 
neglect the effects of neutrals since at T - IQ'^K the chromo- 
spheric plasma is nearly full ionized. Future LRTTM modeling 
with more realistic treatment of spatial structure should include 
the effect of neutrals in deeper cooler regions such as modify- 
ing the resistivity and the collision rate. Collisions were treated 
using a modified form of E with the form i'/v oc v"^, valid for 
V » v,/„ multiplied by the factor v^/(v + v,/,)-' to avoid incorrect 
divergence near thermal speeds v ^ v,/,. 

In the corona, after undergoing nume rous CS £-field accel- 
erations and decelerations, and collisions dTurkmani et alJlToOSL 
[2006) many coronal electrons escaped to the chromosphere, 
a situation geometrically similar to the injection assumed in 
the CTTM with little collisonal HXR emission in the tenuous 
corona. However, once in the chromosphere, as well as colli- 
sions, electrons now undergo re-acceleration by the CS £ fields 
there. This greatly extends some of their lifetimes beyond tcoii, 
increasing the mean photon yield ( and reducing the replenish- 
ment rate T'l hence the beam density ;Fi/Avi needed over area 
A to provide the HXR output J, and so alleviating the problem 
of their large values in the CTTM. In practice the corona ac- 
celerates and injects electrons with a spe ctrum of £, - roughly 
a dou ble power law distribution function dTurkmani et al.ll2005l 
|2006 |). Simulations of these spectral characteristics of injected 
electrons and the resulting bremsstrahlung spectra arising from 
the complex distribution of trajectories E{Et, t) in the thick tar- 
get with re-acceleration (LRTTM)will be the subject of future 
work. Here, to be able to compare simply the dynamics and pho- 
ton yields of electrons arriving in the chromosphere for CTTM 
and LRTTM cases, we limit our analysis to an ensemble of elec- 
trons all of the same 

We have carried out such simulations of E{t) for 10^ elec- 
trons injected randomly with = 50 keV in a chromospheric 
CSC plasma with £, B fields from the MHD simulations dis- 
cussed in Section 3.1 and also for the purely collisonal case. In 
some simulations, chromospheric CS fi values were high enough 
for some electrons in the tail of the local thermal distribution to 



be accelerated but we do not consider these further here. The 
simulation results depend on the electric field which can vary 
from one snapshot to another according to the dynamics of the 
turbulent loop. The main features of our simulations are as fol- 
lows : 

L The accelerating electric field: In the context of particle ac- 
celeration, electric fields are often compared with the Dreicer 
Field &£, (required for the force e&o to overcome collisions 
for a thermal electron of E ^ kT). Fast electrons arriving 
from the corona already have E » kT in the chromosphere 
and the field required for re-acceleration to overcome colli- 
sions for them is smaller than &o by a factor kT/E. 

2. The values of electric fields: The electric fields are zero out- 
side the current sheets and found to take random values in- 
side them. The average of this value in the illustrative case 
used in this paper is £ = 8.2 x 10""* statvolt/cm in the chro- 
mosphere and its maximum value is &max - 2 x 10"^ stat- 
volt/cm. The thickness of the current sheets vary between a 
minimum of 0.5 km and a maximum of 12.5 km. However, 
comparison of cases in terms of 'average' £ values is not 
very meaningful. One could for example have two cases with 
the same volume-averaged £ in one of which £ nowhere ap- 
proached £23 while in the other £ exceeded ££> in some local 
CSs. Differences in the values of £ affect the fraction of the 
injected electrons undergoing re-acceleration, before being 
lost by escape or collisions. They also affect the maximum 
energies electrons reach and their lifetimes. 

3. Direction of the fields: Since the electrons encounter CS £ 
fields in quasi-random directions, the electrons move back 
and forth on guide B fields and their acceleration is stochas- 
tic with E undergoing many changes of sign (cf. Section 2.3) 
as shown in Figure 2a. Though the electric fields are often 
high enough to accelerate or decelerate the electrons inside 
the CSs, there is no global runaway because of the short du- 
rations and quasi-random signs of these kicks. Scattering of 
the electrons also helps enhance their lifetimes by keeping 
them in the CSC region. 

4. Numerical Resistivity and Resolution: The electric field con- 
sidered here is the resistive electric field parallel compo- 
nent and its value depends on the resistivity. The average 
numerical resistivity used here inside a chromospheric CS 
was taken to be 77 = 3 x lO^'-'i (roughly the Spitzer value) 
and, when combined with the numerical resolution used in 
our simulation, results in re-acceleration of electrons in the 
keV-MeV range most relevant to HXR burst production. 
Higher (anomalous) resistivities enhance the re-acceleration 
process. Increasing the resolution of the numerical 3D MHD 
experiment leads to more and thinner fragmented CSs. This 
enhances the re-acceleration process since the electrons un- 
dergo a higher number of smaller kicks, and pass more often 
through electric field free zones in between. 

We found that, in the chromosphere, among the 10^^ test elec- 
tron CSC cases we ran, about 65 % of injected 50 keV electrons 
underwent varying amounts of re-acceleration and lifetime en- 
hancement well beyond the collision time tcoii - 3.2 x lO"-' 
s. Some examples of these re-accelerated electron trajectories 
E(Et, t) are shown in Figure 2a for nine chromospherically re- 
accelerated electrons and for a CTTM electron. It can be seen 
that electrons initially gain and lose energy in the CSs they ran- 
domly pass through, eventually entering electric field free zones 
where they escape or lose their energy to collisions. For our pa- 
rameters, the increase in lifetimes of the injected electrons over 
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tcou ranged from factors slightly higher than unity to around 20 
X with an average over all electrons of about 5 x. 

3.3. Numerical Results for Photon Yield ( in the CSC LRTTM 

As noted in Section 2.3, in the LRTTM scenario, electron trajec- 
tories £■(£■,, t) starting from energy E, are not deterministic but 
stochastic. Thus the only way to arrive at a measure of the pho- 
ton yield f (e, £*) for a single test particle is to use its individual 
equation of motion to compute E(t) for use in time integration 
This is repeated for each of a sample of electrons initially of 
the same energy £, but randomly located then undergoing ran- 
dom kicks in the stochastic £ fields. We want to compare the 
mean ^(e, £*) of these with the CTTM value ^cttmCe, £*)■ In 
the CTTM an electron injected with < e yields no photons 
of energy > e since Ecu < 0. This is not true in the presence of 
re-acceleration since E,„a^- can exceed £». Some criterion there- 
fore has to be adopted for comparison of photon yields. Here 
we chose conservatively to compare the average photon yields 
^(e, £■») for the CTTM and CSC for electrons launched inside 
the source from a specified initial high energy > e namely 50 
keV. 

Based on Equation (|6]l we calcul ated fo r each test electron, 
using the accurate form of Qb from iHaiJgl (1 19971) the rate ^ of 
emission of photons per unit e at 20 keV as a function of t for 
each electron and also the cumulative number emitted up till t 
per unit e as a function of f, hence the total (t —> oo) ^ value 
for each electron - see Figure 2 b, c. In the unique CTTM case 
(Equation ©) for E, = 50 keV CcTTuie, £.) * 2.2 x lO""* pho- 
tons per keV at e = 20 keV. The resulting mean increase in total 
photon yield per keV at 20 keV reached as high as a factor of 
around 20 with an average value around 10. These factors are 
higher than the increases in electron lifetimes mentioned above 
for the reason discussed in Section 2.3. For the given the re- 
sulting increase in yield ^csci^i E*) compared to the collisional 
case increases with e. As we allow e = 10, 20, 30, 40, 45, 49 keV 
to approach £, the relative increase factors Ccsc I Ccttm in aver- 
age yield were about 5, 10, 15, 50, 100, 800. This is because, in 
contrast with the monotonic CCTM fall of E(t) from £», in the 
LRTTM many electrons of initial E - E^ spend many times tc„u 
at £ » e.g. 50 - 800 keV in the example shown in Figure 2. 
This produces many times more photons not only at 20 keV, as 
discussed above, but also at much higher e whereas in the CTTM 
there is no photon yield above 50 keV. A proper comparison of 
LRTTM yield with CTTM is thus rather complicated and will re- 
quire numerical simulations for many and e values. But it is 
clear that the factors quoted above for the LRTTM photon yield 
enhancement are conservative lower limits. 



4. Discussion and Conclusions 

We have shown that substantial re-acceleration in the chromo- 
sphere of electrons accelerated in and injected from the corona 
can greatly reduce the density and number of fast electrons 
needed to produce a HXR burst, and how this might occur in 
a CSC as one example. In the LRTTM. as in the CTTM, most 
electron collisions are in the chromosphere so the LRTTM also 
predicts HXR footpoints. Some of its other properties are, how- 
ever, quite distinct and need much more quantitative work be- 
yond our outline ideas above for the model to be evaluated and 
tested. Here we conclude by briefly discussing some of the issues 
to be addressed. 

1 . Fast Electron anisotropy 



In our CSC simulations we find that the electrons move more 
or less equally up and down the loop axis (< v,+ >^< 
vz- >) with < > / < V|| > about 0.05 in the 
chromosphere and 0.20 in the corona. Unlike the strong 
downward bea ming (< v,+ >»< Vj_ >) in the basic 
CTTM (Brow n.ll972 ). this d istribution is broadly consistent 
with (Kontar & Brow n '2006*) albedo mirror diagnostic 'near 
isotropy' results fromRHESSI spectra. The Vx/vy « 1 prop- 
erty of electrons in the CSC LRTTM may, however, still 
yield enough Hg impact pol arization to contribute to tha t 
observed ( He noux & Chambe 1990; K asparova et alj|2005l) 
though other mechanisms (e.g. fast proton impacts) may also 
contribute (?). 

HXR fine time structure and footpoint synchronism 
When fast electrons in the chromospheric HXR source orig- 
inate by injection from the corona, the HXR light curve 
should reflect the coronal supply rate quite closely since 
even the LRTTM extended fast electron lifetimes t are short. 
So this scen ario is consistent wit h HXR fine time struc- 
ture (< 1) s (Kiplinger et al.' n~983h . footpoint synchronism 
findings (Sakao et al. 1996), and energy -dependent time-of- 
flight delay results OAschwandenl 12004 ). provided that ac- 
celeration in the coronal CSCs is coherent on short enough 
timescales. This coherence should be on the coronal loop 

1/2 

Alfven timescale ta — L/va - O.SL^n^^^ /Bj, so B ^ 500 
gauss suffices to make < Is. 

Interplanetary and HXR Flare Electron Fluxes and Spectra 
In the CTTM the power law spectral index y of HXR emis- 
sion 7(e) is related to the spectral index Sihick of the electron 
injection rate TiEo) by y = 5,iuck - 1 and to the mean source 
electron flux F{E) index 5,,„„ by y = 6,hm + 1 ■ (6thick-6,hin = 2 
because the collisonal energy loss cross section varies as 
E'^). In the LRTTM situation trajectories E{t) are stochas- 
tic and average behaviour depends on the specific CSC re- 
alisation so no such obvious simple relationship exists. This 
complication also means that while integral deconvolution of 
7(e) (e.g. lBrownlll971l; ?) to find the HXR source nVF{E) is 
still fully valid, inference of Tt{E*) is much more difficult 
because of the stochastic character of the electron transport, 
in contrast with the simple CTTM collisional case. 
Using RHESSI and WIND data, ( Krucker et aLll2007l) stud- 
ied the relationship of electron spectra and numbers at the 
Sun to those near the Earth above 50 keV. They find the 
indices and numbers to be well correlated in all events in- 
volving free streaming from the Sun but that the relation- 
ship of spectral indices does not match the CTTM prediction 
(Krucker et al. 2007, 2009). Further, the numbers of elec- 
trons in IP space an d in Type III Burst s are smaller by a 
factor of order 500 dKrucker et al.l 120071) than required for 
the HXR source in the CTTM model. The numbers required 
for microwave bursts are also often found to be considerably 
less than for HXRs in the CTTM interpretatio n though this 
number is very sensitive to assumed conditions dLee & Gary! 
(2000) . For these to be consistent with the LRTTM reduc- 
tion in electron numbers would imply even more effective 
re-acceleration than in the illustrative example we gave here, 
such as due to anomalous resistivity. 
Impulsive flare heating 

Various aspects of impulsive flare heating data have been in- 
voked in support of the CTTM, including the Neupert Effect 
that flare soft XR li ght curves cor relate with the integral of 
HXR light curves (lNeupertlll968l) . This is often attributed 
to CTTM collisional heating of the SXR source plasma by 
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HXR emitting fast electrons. However, the observed relative 
time sequences of EM{t), T(t) are hard to reconcile with this 
in any obvi ous way even when filamented loop struct ures are 
considered dVeronig et al.ll2005l: [Stoiser et al.ll2008bl) . In the 
LRTTM the total power delivered to fast electrons in order 
to offset collisional losses is comparable to that in the CTTM 
model and can likewise heat the impulsive flare atmosphere, 
though the spatial distribution of that heating can be very 
different from the CTTM case. In the proposed LRTTM sce- 
nario the coronally injected beam rate fi is reduced consid- 
erably from the CTTM rate so beam heating of the corona 
is reduced from its CTTM value. However, for given HXR 
output, the total beam power involved in the whole HXR 
source has to be at least as large as in the CTMM. Thus in 
the LRTTM more power goes into chromospheric heating as 
re-acceleration drives fast electrons against collisional losses 
there. In addition, if the extended electron lifetimes result in 
their penetrating deeper, beam heating may be effective to 
much greater depths than in the CTTM. This might offer a 
solution to the problem of heating white light flares by elec- 
tron beams (Neidig 1989; Fletcher et al. 2007). 
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